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Abstract 

A  method  to  recover  atmospheric  temperature  profiles  using  a  ground-based 
Fourier  Transform  Infrared  spectrometer  was  investigated.  The  method  used  a  difference 
form  of  the  radiative  transfer  equation,  a  Bomem  MR  series  Fourier  Transform  Infrared 
Spectrometer  to  collect  atmospheric  radiance  values,  and  the  Phillips  Laboratory  Expert- 
assisted  User  Software  (PLEXUS)  atmospheric  radiance  model,  to  recover  an 
atmospheric  temperature  profile.  The  method  researched  uses  radiance  values  from  both 
the  spectrometer  measurements  and  the  atmospheric  model,  along  with  kernel  functions 
calculated  by  the  atmospheric  model  as  input  to  a  difference  form  of  the  radiation  transfer 
equation.  From  this  the  change  in  brightness  temperatures  was  determined.  The  method 
assumes  that  the  actual  brightness  temperature  profile  is  a  summation  of  a  standard  or 
reference  brightness  temperature  profile  plus  some  change  in  the  brightness  temperature. 
The  brightness  temperature  profile  used  by  the  atmospheric  model  is  the  reference 
brightness  temperature  profile.  Planck’s  Law  was  employed  to  transform  the  calculated 
brightness  temperature  function  into  a  temperature  function.  A  temperature  profile  was 
retrieved,  although  significant  differences  existed  between  the  recovered  temperature 
profile  and  a  radiosonde  recovered  temperature  profile. 


XI 


ATMOSPHERIC  TEMPERATURE  PROFILES  BY  GROUND-BASED  INFRARED 


SPECTROMETER  MEASUREMENT 
I.  Introduction 

1.1  Summary 

In  the  last  decade,  great  strides  have  been  made  in  the  area  of  remote  sensing  by 
surface  based  equipment.  Advances  in  technology  and  reduction  in  the  cost  associated  in 
the  manufacture  of  remote  sensing  equipment  have  made  the  plausibility  of  remote 
sensing  a  reality.  In  particular,  remote  sensors  are  sought  for  as  a  means  of  replacing  or 
accenting  current  atmospheric  measurement  techniques. 

Radiosondes  are  the  principle  devices  used  to  obtain  atmospheric  profiles.  A 
single  radiosonde  costs  approximately  $100.  Typically,  radiosonde  sites  launch 
radiosondes  twice  a  day.  Even  with  these  conservative  costs,  a  yearly  supply  of 
radiosondes  for  one  station  can  cost  over  $73,000.  This  estimate  does  not  include  the  cost 
of  tracking  the  radiosonde  or  disseminating  the  data.  Using  surface-based  remote  sensors, 
atmospheric  profiles  can  be  recovered  (Smith  et  al.  1990).  The  infrared  spectrometer  is 
one  instrument  being  studied  for  its  use  in  obtaining  infrared  spectral  data  to  produce 
atmospheric  profiles.  Smith  et  al.  (1990)  believes  a  basic  interferometer  system  capable 
of  performing  ground-based  atmospheric  profiling  would  cost  less  than  $50,000.  A 
ground-based  infrared  spectrometer  system  offers  the  possibility  of  saving  millions  of 
dollars  annually. 

A  problem  with  the  current  profiling  system,  radiosondes,  is  the  temporal 
resolution  of  upper  air  profiles.  Typically,  radiosondes  are  launched  twice  a  day,  at  OOZ 
and  12Z,  and  recover  one  atmospheric  profile  for  each  radiosonde  launched.  Because 
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atmospheric  profiles  of  the  troposphere  are  constantly  changing,  the  recovery  of 
atmospheric  profiles  twice  a  day  leaves  large  temporal  gaps  in  the  data.  An  infrared 
spectrometer  has  the  ability  to  produce  atmospheric  radiance  measurements  every  few 
seconds.  From  these  radiance  measurements,  a  profile  can  be  calculated  for  each  radiance 
measurement.  An  increased  temporal  resolution  provided  by  ground-based  spectrometers 
will  enhance  all  areas  of  meteorology. 

With  the  cost  and  the  temporal  coverage  benefits,  the  use  of  an  infrared 
spectrometer  to  supplement  or  replace  the  radiosondes  as  the  primary  mean  of  retrieving 
atmospheric  profiles  looks  promising.  The  remaining  question  is  which  process  to  use  to 
recover  atmospheric  profiles. 

1 .2  Statement  of  Problem 

This  research  used  a  difference  form  of  the  radiative  transfer  equation,  a  Bomem 
MR  series  Fourier  Transform  Infrared  Spectrometer  (FTIR)  to  collect  atmospheric 
radiance  values,  and  the  Phillips  Laboratory  Expert-assisted  User  Software  (PLEXUS) 
atmospheric  radiance  model,  to  recover  an  atmospheric  temperature  profile.  The  method 
researched  uses  radiance  values  from  spectrometer  measurements  and  the  atmospheric 
model,  along  with  transmissivity  value  calculated  by  the  atmospheric  model  as  input  to 
the  difference  form  of  the  radiation  transfer  equation.  From  this  the  change  in  brightness 
temperatures  can  be  determined.  The  method  assumes  that  the  actual  brightness 
temperature  profile  is  a  summation  of  a  standard  or  reference  brightness  temperature 
profile  plus  some  change  in  the  brightness  temperature.  The  brightness  temperature 
profile  used  by  the  atmospheric  model  is  the  reference  brightness  temperature  profile. 
The  goal  is  to  determine  the  changes  in  the  atmospheric  brightness  temperature  profile 
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and  to  recover  the  actual  brightness  temperature  profile.  The  temperature  profile  can  be 
calculated  using  the  relationship  between  brightness  temperature  and  temperature  in 
Planck’s  Law. 

1 .3  Overview 

This  thesis  consists  of  five  chapters.  Chapter  II  contains  literature  review 
material,  which  contains  fundamental  material  to  support  this  research.  Chapter  III 
describes  the  procedures  used  during  this  research.  Chapter  IV  presents  the  results  of  this 
research.  Chapter  V  discuses  conclusions  made  from  this  research. 
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II. 


Literature  Review 


2.1  Atmospheric  Radiation  Theory 

The  concept  of  using  radiative  transfer  to  formulate  the  absorbing,  scattering,  and 
creation  of  natural  radiation  within  a  volume  filled  with  particles  interacting  with 
radiation  was  developed  for  astrophysical  problems  (Schanda,  1986).  Radiation  can  be 
absorbed  by  atmospheric  gases  in  continuum:  (ionization  and  photo-dissociation),  and 
discrete  (electronic  transitions,  vibrational  transitions,  and  rotational  transitions)  bands. 

In  the  infrared  region,  vibrational  and  rotational  transitions  account  for  most  of  spectral 
line  emissions  (Goody  and  Young,  1989).  Transitions  between  two  energy  levels  produce 
discrete  spectral  line  emission  or  absorption,  but  due  to  molecule  collision,  spectral 
broadening  occurs  and  distorts  a  molecule’s  energy  levels.  This  distortion  in  the  energy 
level  causes  the  transitions  to  occur  over  a  range  of  frequencies  instead  of  one  discrete 
frequency. 

The  radiative  transfer  equation  (RTE)  accounts  for  all  four  possible  radiative 
interactions:  the  emission  of  radiation  by  a  gas,  the  absorption  of  radiation  by  a  gas,  the 
scattering  of  radiation  out  of  a  beam,  and  the  scattering  of  radiation  into  a  beam.  A 
general  formula  for  the  instantaneous  change  of  the  radiance  is: 

—  =A+B+C+D  (2.1) 

& 

where, 

L  =  monochromatic  radiance  (W/cm  cm'  sr) 

s=  distance  along  a  beam  (meters) 
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A  =  absorption  of  radiation  by  a  gas  (W/cm^  cm"’  m"') 

B  =  emission  of  radiation  by  a  gas  (W/cm^  cm"’  m"’) 

C  =  scattering  of  radiation  out  of  a  beam  (W/cm^  cm"’  m"’) 

D  =  scattering  of  radiation  into  a  beam  (W/cm^  cm"’  m"’) 

A  special  case  of  the  RTE  is  sometimes  referred  to  as  the  Schwarzchild's 
equation.  The  Schwarzchild  equation  assumes  no  scattering  of  radiation  into  or  out  of  the 
beam.  On  a  clear,  cloudless  day  the  Schwarzchild's  equation  (Wallace,  1997),  Equation 

2.2,  is  valid  for  the  atmosphere,  and  is: 

00 

L,=I^+\b,{z)W,{z,/2^z  (2.2) 

0 

where, 

9  1 

=  spectral  radiance  at  the  top  of  the  atmosphere  (W/cm  cm"  sr) 

9  1 

=  spectral  radiance  received  (W/cm  cm"  sr) 
ju  =cos(0) 

0  =  zenith  angle 

(z)  =  blackbody  radiance  (W/cm  cm" ) 

{z,/j)  =  kernel  function  (change  of  transmittance  with  respect  to  height) 

z  =  height  (meters) 

X  =  wavelength  (meters) 

Retrieval  approaches  use  an  inverse  form  of  the  radiative  transfer  equation  to 
solve  for  atmospheric  temperature  profiles.  Using  the  Schwarzchild's  form  of  the 
radiative  transfer.  Equation  (2.2),  and  calibrated  radiance  values  obtained  from  an 
infrared  spectrometer,  it  is  possible  to  calculate  a  temperature  profile.  This  version  of  the 
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radiative  transfer  equation  assumes  no  scattered  or  diffuse  radiation.  These  assumptions 
are  viable  for  the  infrared  spectrum  in  the  atmosphere  in  the  absence  of  clouds  (Kidder 
and  Vonder  Haar,  1995). 

A  linear  approximation  for  the  radiative  transfer  equation  in  a  matrix  form  is 
(Schanda,  1986): 

=  tf'.B,  (2.3) 

where, 

=  vector  of  monochromatic  radiance  values  (W/cm^  cm’'  sr) 

=  a  matrix  of  change  in  transmittance  values 

=  vector  of  brightness  temperatures  (W/cm  cm'  ) 

Using  the  ground-based  spectrometer,  the  measured  radiance,  Lx,  is  observed  and 
the  weighting  functions,  W^(z,//) ,  can  be  calculated  using  an  atmospheric  model.  An 

inverse  method  uses  the  observed  radiances  to  calculate  a  Planck's  radiance  profile.  The 
inverse  of  the  non-scattering  radiative  transfer  equation  is: 

IV;'Z,=B,  (2.4) 

With  Bx  known,  then  temperature  is  calculated  using  Planck’s  function. 

=C,A-'  /[expiC^/ATJ-l] 

(2.5) 

where, 

2  1 

Bx  =  brightness  temperature  (W/cm  cm' ) 
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T  =  temperature  (Kelvin) 

X  =  wavelength  (meters) 

Cl  =  1.1919439  X  10''^  Wm^ 

C2  =  1.438769  X  10'^  mK 
2.2  Retrieval  Approaches 

In  late  1988,  researchers  conducted  the  Ground-based  Atmospheric  Profiling 
Experiment  (GAPEX)  (Smith  et  al,  1990).  One  part  of  GAPEX  involved  using  a  high- 
resolution  interferometer  (HIS)  to  measure  the  infrared  spectrum.  Atmospheric 
temperature  and  water  vapor  profiles  were  calculated  using  the  measured  radiances  from 
the  HIS.  Researchers  used  the  inverse  of  the  non-scattering  radiative  transfer  matrix, 
Equation  2.4,  to  recover  temperature  profiles.  The  HIS  was  used  to  obtain  atmospheric 
radiance  data  and  statistical  data  was  used  to  determine  the  kernel  function. 

In  mid  1992,  researchers  experimented  with  retrieving  atmospheric  profiles  firom 
radiance  values  obtained  by  a  Double  Beam  Interferometer  Sounder  (DBIS)  (Theriault 
1993);  (Theriault  et  al,  1996).  The  DBIS  experiment  used  the  approach  based  on  the 
minimum  information  method.  There  are  two  main  advantages  to  this  approach:  one, 
"first-guess"  profiles  are  not  required,  and  second,  simultaneous  retrieval  for  multiple 
parameters  is  possible. 

The  method  works  as  follows:  first,  researchers  assumed  that  the  measured 
radiance  is  a  function  of  the  actual  radiance  plus  a  linear  perturbation  (Theriault  1993); 
(Theriault  et  al,  1996).  The  perturbations  are  linearized  about  a  point  Lq, 

L(v)  =  L,{v)  +  dL(v)  (2.6) 
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where  L(v)  is  the  measured  spectral  radiance,  Lo(v)  is  the  actual  radiance,  and  dL(v)  is 
the  perturbation  of  the  radiance.  Second,  the  approach  employs  a  linear  difference 
technique  in  a  matrix  form. 

B,^[w^s:^w+si'rw^s:'L,  (2.7) 

This  approach  constrains  the  solution  by  weighting  the  equation  with  two  error  matrices, 
Se  and  Sb-  The  error  weighting  matrices  accounts  for  instrument  error  and  uncertainty  in 
the  "first-guess"  profile. 

Obtaining  the  damping  matrix,  Sb,  poses  a  problem  (Theriault  1993:);  (Theriault 
et  al,  1996).  The  selection  of  the  damping  matrix  must  be  completed  empirically. 
Another  problem  deals  with  the  time  required  to  iterate  to  a  solution. 

In  another  experiment,  researchers  at  the  University  of  Wisconsin-Madison  use 
the  Atmospheric  Emitted  Radiance  Interferometer  (AERJ)  to  measure  atmospheric 
radiance  every  ten  minutes.  Atmospheric  temperature  profiles  were  produced  using  the 
measured  radiance.  AERI  uses  an  inverse  technique  to  retrieve  a  temperature  profile 
from  ground-based  radiance  measurements.  Using  the  inverse  of  the  radiative  transfer 
equation,  researchers  at  the  University  of  Wisconsin-Madison  use  a  statistical  first  guess 
of  an  atmospheric  temperature  profile  and  the  measured  radiances  to  solve  the  inverse 
radiative  transfer  equation.  The  University  of  Wisconsin-Madison  has  detailed 
information  on  the  AERI  project  at  http://cimss.ssec.wisc.edu/aerlwww/aeri/. 

Hauser  (1999)  experimented  with  the  Bomem  MR- 154  FTIR  spectrometer  to  recover 
infrared  radiance  measurements  of  the  atmosphere.  The  retrieval  method  used  a  least 
squares  solution  to  obtain  a  temperature  profiles  from  observed  radiances.  The  method 
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involves  multiplying  both  sides  of  equation  2.2  by  the  transpose  of  the  weighting 
function,  which  Hauser  identifies  as  the  “kernel”,  and  solving  for  blackbody 
temperatures.  The  weighting  functions  or  kernels,  W,  are  the  changes  in  monochromatic 
transmissivity  over  path  length.  MODTRAN  calculations  are  used  to  compute  the 
weighting  functions.  Once  the  weighting  function  was  determined,  the  inverse  matrix 
was  solved. 

=  B,  (2.8) 

In  his  results,  Hauser  (1999)  states  that  the  procedure  did  not  work,  but  gives  no 
explanation  as  to  why. 

2.3  Atmospheric  Models 

2.3.1  PLEXUS 

PLEXUS  is  a  user-oriented  software  platform,  which  integrates  user  input  into 
DOD  standard  atmospheric  and  celestial  models.  PLEXUS  provides  a  graphical  interface 
for  the  following  Air  Force  Research  Lab  (AFRL)  codes:  MODTRAN3.7,  SAMMl, 
FASCODE3P,  and  SAGl .  The  PLEXUS  software  allows  the  user  to  input  parameters, 
then  the  PLEXUS  software  converts  the  numeric  input  to  the  native  units  of  the 
applicable  AFRL  code.  The  source  codes  of  the  models  are  not  altered  in  any  way. 
PLEXUS  is  available  on  a  CD  through  AFRL  (PLEXUS,  1999). 
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Spectral  Radiance  (W/cnn2-sr-cm-1 ) 


Frequency  (cm-1) 


Figure  2.1  Example  of  graphical  radiance  output  from  PLEXUS. 


Frequency  (cm-1) 


Figure  2.2  Example  of  graphical  transmittance  output  from  PLEXUS. 
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2.3.2  SHARC  and  MODTRAN  Merged  (SAMM) 


The  SAMM  model  was  used  to  calculate  all  model  derived  transmittance  and 
radiance  data.  SAMM  is  an  atmospheric  model,  which  combines  the  Moderate 
Resolution  Transmittance  Model  (MODTRAN)  and  the  Strategic  High- Altitude 
Atmospheric  Radiance  Code  (SHARC).  MODTRAN  was  developed  by  the  Air  Force 
Research  Laboratory  as  a  way  to  calculate  path  radiance  and  path  transmittance  in  the 
infrared,  visible,  and  near  ultraviolet  spectral  regions  for  a  given  atmospheric  path  below 
100  km.  MODTRAN  is  a  band  model  that  uses  empirically  derived  data  to  compute 
transmittance  and  radiance  values.  SHARC3  is  a  non-Local  Thermodynamic  Equilibrium 
(NLTE)  code  for  computing  infrared  path  radiance  for  arbitrary  paths  between  50  and 
300  km  in  the  2-40  micron  spectral  region.  It  also  incorporates  SHARC  Atmospheric 
Generator  (SAG).  SAG  is  an  AFRL  code  that  calculates  atmospheric  profiles  based  on 
geophysical  and  geographical  information.  SAG  uses  the  NRL  climatological  database 
(PLEXUS,  1999),  which  is  comprised  of  atmospheric  profiles  of  mean  monthly 
concentrations  at  1  to  5  km  increments  and  ten-degree  latitude  increments.  SAG 
interpolates  between  these  values  and  PLEXUS  uses  the  output  for  SAMM. 

2.4  Bomem  MR- 154  Spectrometer 

The  Bomem  MR- 154  Fourier  Transform  Infrared  (FTIR)  Spectoradiometer  is  an 
instrument  capable  of  obtaining  highly  accurate  spectral  features  of  a  target  source  in  the 
525  cm’’  to  6000  cm"'  (19  pm  to  5.55  pm)  range  with  a  maximum  resolution  of  1  cm'^ 
The  MR- 154  FTIR  spectrometer  is  comprised  of  two  subsystems:  the  data  processing 
and  control  assembly  and  a  Michelson  interferometer. 
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The  data  processing  and  control  system  consists  of  a  personal  computer,  which 
runs  the  Research  Acquire  software  package.  The  software  package  provides  a  user 
interface  to  control  the  interferometer  to  collect  data.  The  software  also  allows  a  user  to 
process  the  data  to  produce  calibrated  radiance  files. 

Figure  2.3  is  a  top  view  of  the  interferometer  with  the  main  components 
identified.  The  main  components  consist  of  two  detectors,  a  collimator  with  an  aperture, 
transfer  optics  with  apertures  and  a  cold  subtraction  source.  The  interferometer’s  two 
detectors  are  the  mercury,  cadmium,  and  telluride,  MCT,  a  solid-state  detector  and  the 
indium  antimonide,  InSb,  a  solid-state  detector.  The  MCT  detector  allows  for  the 
collection  of  radiometric  data  in  the  525  cm’’  to  1800  cm*  (19  pm  to  5.55  pm)  region. 
The  InSb  detector  allows  for  the  collection  of  radiometic  data  in  the  1 800  cm  *  to  6000 
cm'*  (5.56  pm  to  1.67  pm)  region.  The  cold  subtraction  source  and  the  detectors  are  filled 
with  liquid  nitrogen  to  reduce  the  thermal  noise  of  the  instrument.  The  collimator  and 
transfer  optics,  along  with  their  apertures,  control  the  radiometric  data  reaching  the 
detectors. 

Figure  2.4  is  reproduced  from  Hauser  (1999).  It  displays  the  principle  of  operation 
of  the  FTIR  spectrometer.  Source  radiation  enters  the  spectrometer  through  the  input 
collimator  (1).  Then  radiation  is  split  by  a  beam  splitter,  with  one  beam  reflected  by  a 
stationary  mirror  while,  the  other's  path  length  is  varied  with  a  moving  mirror  (Dl).  Then 
the  radiation  beam  is  recombined  (2)  and  recorded  as  an  interferogram.(3)  A  Fourier 
transform  is  applied  to  the  interferogram  to  produce  a  raw  spectrum  (4).  The  raw 
spectrum  is  calibrated  using  a  previously  determined  calibration  spectrum  (5).  The  final 
output  is  calibrated  spectral  data  (6).  Further  information  on  the  components  of  the 
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Bomem  MR- 154  FTIR  Spectoradiometer  and  its  operation  can  be  found  in  the  Bomem 
MR  Series  Documentation  set  (BOMEM  Inc.  1997). 


Transfer  optics  and  aperture  Cold  subtraction  source 


Figure  2.3  Bomen  MR-154  Basic  Components. 
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Top  View  of  IntttTial  Components 


Source 


1 


3 


6 


Recombined 


Fourier 

Transform 

Applied 
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Calibrated  Spectral 


► 
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Calibration 


Figure  2.4.  Principle  of  operation  of  the  FTIR  spectrometer.  Incident  energy,  1, 
aperture  and  encounters  beam  splitter  at  point  2.  The  beam  is  recombined  where 
and  destructive  interference  result  represented  by  the  interferogram  at  3.  A  Fourier 
is  applied  to  the  interferogram  and  a  raw  spectrum  results.  At  5,  a  calibration 
completed.  Finally,  a  calibrated  spectral  plot  is  obtained  at  point  6  (after  BOMEM  Users 
Manual.)  (Hauser,  1999). 
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III.  Methodology 


3.1  Overview 

This  chapter  presents  the  methodology  used  during  the  research.  Starting  with 
Equation  2.2,  Schwarzchild's  equation  for  the  actual  radiance  is: 

L,^\w,{Z)B,{z)dz  (3.1) 

where, 

=  Monochromatic  radiation  for  a  specific  wavelength  of  the  atmosphere 
(W/cm^  cm''  sr'') 

ff((z)  =  Kernel  function  for  a  specific  wavelength  of  the  atmosphere 

i^(z)  =  Brightness  temperature  function  for  a  specific  wavelength  of  the 
atmosphere  (W/cm^  cm  ') 

z  =  height  (meters) 

A  second  radiation  equation  is  employed  as  a  reference: 

I,  =  jmZ)B,(z)dz  (3.2) 

where, 

Zj  =  Monochromatic  radiation  for  a  specific  wavelength  from  an  atmospheric 
model  (W/cm^  cm''  sr'') 

ff^(z)  =  Kernel  function  for  a  specific  wavelength  from  an  atmospheric  model 

B^(z)  =  Brightness  temperature  function  for  a  specific  wavelength  from  an 
atmospheric  model  (W/cm  cm' ) 

Subtracting  Equation  3.2  from  Equation  3.1  provides  a  relationship  between  the 
atmospheric  model  derived  parameters  and  the  atmospheric  parameters. 
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Z,  -4  =  I  Wiz\mz)dz-  J  W{z\BJiz)dz 


(3.3) 


Assuming 


^(z)=Z^(z)+A5(z)  (3.4) 

using  A6(z)  =Change  in  brightness  temperature  function  (W/cm^  cm’’)  then: 

z,  -4 = j  w{z\{B,{^)+mz))dz-  j  mzrB,{z)dz  (3.5) 

Expanding  terms  and  rearranging, 

4  -4  =  J  mz)B,{z)dz+  J  W{z\mz)dz-  J  W{z\  *B,{z)dz  (3.6) 

4  -4  -  J  W,(z)B,(z)dz+ 1  W(z\  *B,{z)dz=  J  W{z\mz)dz  (3.7) 


If  the  kernel  functions,  W(z),  that  do  not  change  with  time  can  be  found,  that  is 
Wi(z)  =  W2(z)  then. 


4-4  =  JlF(z)AS(z)d'z 


(3.8) 


In  this  research,  W,  L2  and  B2  are  determined  using  PLEXUS  and  its  integrated 
atmospheric  model,  SAMMl.  Li  is  actual  radiometric  measmements  using  the  Bomem 
MR- 154  spectrometer. 

We  construct  an  algorithm  to  find  the  temperature  from  the  radiometric 
measurements. 

Step  1.  We  first  search  regions  of  the  spectrum  where  the  kernel  functions  are  insensitive 
to  the  state  of  the  atmosphere  while  there  is  a  significant  change  in  the  observed  radiance 
(see  Section  3.2). 
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Step  2.  Kernel  functions  are  calculated  using  the  PLEXUS  atmospheric  model.  The 
kernel  functions  will  be  determined  by  calculating  the  transmittance  at  a  number  of 
specified  altitudes  and  then  determining  the  path  derivatives  dx/dz  (see  Section  3.3). 

Ideally,  the  kernel  functions  should  be  continuous,  because  the  change  in 
transmittance  is  continuous  in  the  atmosphere.  The  kernel  functions  calculated  with 
PLEXUS  are  discrete.  For  this  reason,  continuous  functions  are  fitted  to  the  discrete 
kernel  values  (see  Sections  3.3  and  3.4).. 

PLEXUS  uses  predetermined  atmospheric  profiles  to  calculate  radiance  outputs. 
Kernel  functions  were  tested  by  comparing  a  recovered  temperature  profile,  using  the 
PLEXUS  radiance  output  and  previously  determined  kernel  functions,  to  the  model  input 
temperature  profile. 

Step  3.  The  final  step  is  to  use  the  spectrometer  to  recover  actual  radiance  measurements 
of  the  atmosphere  and  determine  the  temperature  profile  using  the  calculated  kernel 
functions.  This  resulting  temperature  profile  will  be  compare  to  a  local  radiosonde- 
recovered  temperature  profile  coincident  in  time  (see  Sections  3.5  and  3.6).. 

3.2  Spectral  Region  Determination 

Spectral  regions  from  500  cm"'  to  3000  cm''  were  chosen  based  on  the  Bomem 
spectrometer’s  spectral  range.  To  minimize  the  number  of  possible  wavelengths,  regions 
of  the  spectrum  were  disqualified  as  being  poor  choices  for  recovering  tetftperature 
profiles.  The  first  was  to  disqualify  all  regions  whose  absorption  is  mainly  due  to  highly 
variable  density  profiles  of  atmospheric  constituents.  Goody  and  Yung  (1989)  list  three 
constituents  as  variable  or  highly  variable  and  give  the  appropriate  absorption  bands 
associated  with  the  constituents.  Table  3.1  lists  these  constituents  and  their  absorption 
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bands.  These  absorption  bands  were  disqualified.  Disqualifying  these  absorption  bands 
leaves  the  band  from  2400  cm  ’  to  2800  cm’’. 


Table  3.1  Highly  Variable  Atmospheric  Constituents  and  their  Absorption  Bands. 


Constituent 

Absorption  Band 

Transition 

H20,  water 

0  ~  1000  cm-1 

rotational 

H20,  water 

900  ~  2400  cm-1 

vibrational 

H20,  water 

2800  ~  4400  cm-1 

vibrational 

03,  ozone 

1000-  1100  cm-1 

vibrational 

CO,  carbon  monoxide 

2100 -2400  cm-1 

vibrational 

Two  dates  were  chosen  based  on  the  National  Research  Laboratory,  NRL,  data 
base  temperature  profiles.  These  temperature  profiles  were  used  as  the  inputs  to 
PLEXUS.  The  months  of  January  and  July  were  chosen  based  on  the  greatest  differences 
in  the  temperature  profiles.  PLEXUS  software  was  used  to  calculate  radiance  value  for 
650  cm'*  to  3000  cm"’.  Figure  3.1  is  the  PLEXUS  derived  January  radiance  profile,  while 
Figure  3.2  is  the  PLEXUS  derived  July  radiance  profile.  The  values  are  calculated  for  a 
path  length  of  infinity,  the  total  atmosphere.  Figure  3.3  is  the  ratio  of  the  two  PLEXUS 
derived  radiance  profiles. 

Transmittance  values  were  calculated  using  PLEXUS  over  the  650  cm"’  to  3000 
cm'*  range  for  both  January  and  July.  At  each  discrete  height,  the  quotients  of  the  July 
and  January  transmittance  values  were  determined.  Wavenumbers  where  the  quotients 
were  one  or  near  one  and  there  was  a  significant  change  in  the  radiance  were  considered 
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Spectral  Radiance  (W/cm2'Sr-cm-1) 


to  be  caused  by  mostly  temperature  changes  in  the  input  atmospheric  profiles.  Figure  3.4 
displays  the  ratio  of  the  January  and  July  transmittance  values  calculated  for  a  path  length 
of  the  entire  atmosphere. 


Frequency  [cm-1) 


Figure  3.1  PLEXUS  graphical  radiance  output  for  15  January  2000. 


Figure  3.2  PLEXUS  graphical  radiance  output  for  15  July  1999. 
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3.3  Calculating  Kernel  Functions 

PLEXUS  was  used,  in  a  batch  mode,  see  Appendix  2,  to  calculate  transmittance 
values  over  the  650  cm’’  to  3000  cm  *  range  with  varying  pathlengths.  The  following 
pathlengths  were  used,  see  Table  3.2.. 


Table  3.2  Pathlengths  used  in  calculation  of  transmittance. 


Lower  Height 

Upper  Height 

Change  in  Pathlength 

0  meters 

2  kilometers 

10  meters 

2  kilometers 

10  kilometers 

50  meters 

10  kilometers 

50  kilometers 

1000  meters 

50  Kilometers 

300  kilometers 

5000  meters 

For  selected  wavenumbers,  the  following  equation  form.  Equation  3.9,  was 
calculated  by  a  best- fit  software  package  to  fit  a  function  to  the  discrete  transmittance 
values. 


The  equation  form.  Equation  3.10,  was  used  for  all  calculated  transmittance 
functions,  x(z)  as  a  function  of  height,  z. 


z(z)  = 


fl+CZ 

\+bz+di 


(3.9) 


Next,  the  derivative  of  transmittance  with  respect  to  height,  that  is,  the  kernel 
function,  was  determined.  Equation  3.14  shows  the  derivative  of  x(z)  and  was  used  for  all 
used  for  all  calculated  kernel  functions,  W. 

-{-c+cdi  +atH-2ad^ 


W(z) = di(z)  /  dz— - 


(3.10) 


The  parameters  for  the  transmittance  and  kernel  functions  calculated  can  be  found  in 
Appendix  1. 
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The  following  is  an  example  of  the  method  used  to  determine  the  kernel 


functions.  Table  3.3  is  an  example  of  transmittance  values  at  given  heights.  Figure  3.5  is 
the  output  from  the  best- fit  software  using  data  from  Table  3.3. 

Table  3.3  Example  of  transmittance  values  for  specific  pathlengths. 


Path  Length  (m) 
0 

100 

200 

300 

400 

500 

600 

700 

800 

900 

1000 


Transmittance 

1.000000 

0.816010 

0.714550 

0.642030 

0.584930 

0.538240 

0.499040 

0.465520 

0.436640 

0.411210 

0.389160 


y=(a+cx)/(1+bx+dx2) 

a=0.99832304  b=0.0084152028  c=0.0057232514  cl=7.8675985e-06 


Figure  3.5.  Example  of  output  from  best-fit  software  using  values  from  Table  3.3. 
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3.4  Calculating  Atmospheric  Brightness  Temperature  Functions 

Atmospheric  brightness  temperatures  were  calculated  using  Planck’s  formula  and 
physical  temperature  provided  by  the  atmospheric  model.  Temperatures  from  the 
PLEXUS  software  were  given  at  discrete  height  values,  see  Figure  3.6. 


Figure  3.6  Temperature  profile  used  by  PLEXUS  for  15  November  1999. 

These  temperature  values  were  used  to  calculate  the  brightness  temperature  at 
specific  wavelengths  using  Planck’s  law. 

=  qX'  /[exp(C2  /AT)  - 1]  (3.11) 

where, 

Bx  =  brightness  temperature  (W/cm^  cm*') 

T  =  temperature  (Kelvin) 
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X,  =  wavelength  (meters) 

Ci=  1.1919439  X  10''®  Wm^ 

C2  =  1.438769  X  10'^  m 

Once  the  discrete  brightness  temperature  profile  is  determined,  a  best- fit  software  is  used 
to  fit  a  function  to  the  data.  The  brightness  temperature  function  calculated  for  this 
research  can  be  found  in  Appendix  1 . 

The  following  is  an  example  of  the  method  used  to  determine  the  brightness 
temperature  functions.  Table  3.4  is  an  example  of  calculated  brightness  temperature  at 
given  heights  and  temperatures  using  a  wavelength  of  2568  cm"'.  Figure  3.7  is  the  output 
from  the  best-fit  software  using  data  from  Table  3.4. 

Table  3.4  Example  of  temperature  and  calculated  brightness  temperatures  for  2568  cm'\ 


Height  (m) 

Temperature  (K) 

0 

290.2 

1000 

286.9 

2000 

281.9 

3000 

275.6 

4000 

268.4 

5000 

260.9 

6000 

253.4 

7000 

246.0 

8000 

239.2 

9000 

232.9 

10000 

227.4 

11000 

222.8 

12000 

219.0 

13000 

216.0 

14000 

213.7 

15000 

212.1 

16000 

211.3 

17000 

211.0 

18000 

211.1 

19000 

211.4 

20000 

211.9 

2  1 

Brightness  Temperature  (W/cm  cm"  ) 
3.91E-07 
3.38E-07 
2.68E-07 
1.98E-07 
1.39E-07 
9.34E-08 
6.12E-08 
3.96E-08 
2.57E-08 
1.70E-08 
1.16E-08 
8.27E-09 
6.20E-09 
4.89E-09 
4.07E-09 
3.58E-09 
3.34E-09 
3.26E-09 
3.29E-09 
3.39E-09 
3.53E-09 
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y=  (a +c  X +eS+g  +i  +l<)?)/(1 +b  X +d )? +f)e3 +h  )(^ +j  >?) 
a=3,9076456e-07  b=-0.00020002041  c=-1.1778633e-10  d=5.958575e-08 
■  e=1 .58601 990-1 4 1i:-8.3255991  e-1 2  g=-1 .1 397826e-1 8  h=7.998501 6e-1 6 
i=4.30601 940-23  ^-2.74384580-20  k=-6.7821 1 01  e-28 


Figure  3.7  Best- fit  brightness  temperature  function  using  Table  3.4  data. 

3.5  Collecting  Atmospheric  Radiance  Measurements 

The  MR- 154  spectrometer  was  used  to  obtain  atmospheric  radiance 
measurements.  The  following  sections  details  the  setup  of  the  equipment,  equipment 
calibrations  and  the  recovering  of  atmospheric  radiance  values. 

3.5.1  Spectrometer  Setup 

The  setup  of  the  MR- 154  spectrometer  was  conducted  following  the  procedures  in 
the  Bomem  MR  Series  Documentation  Set  (1997).  The  MR-154  spectrometer’s  was 
deployed  on  the  rooftop  of  Bldg  640  at  Wright-Patterson  AFB.  The  spectrometer's  tripod 
was  set  up  and  leveled  using  the  bubble-level  located  on  the  tripod.  Once  the  tripod  was 
leveled  the  azimuth  was  calibrated  using  a  compass.  Next  the  main  assembly  was 
attached  to  the  tripod  and  the  cold  subtraction  source  and  detectors  were  attached  to  the 
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main  assembly.  The  cold  subtraction  source  and  detectors  were  filled  with  liquid 
nitrogen  to  ensure  a  minimal  thermal  noise.  The  collimator,  amplifier  and  resolution 
settings  were  set  to  the  best  setting  determined  by  Hauser  (1999):  the  gain  was  full,  the 
aperture  was  at  diameter  of  6.4  mm,  and  the  resolution  was  1  cm''.  The  cover  to  the  main 
assembly  was  attached  and  the  level  of  the  main  assembly  was  rechecked.  The  main 
assembly’s  internal  heater  was  set  to  30  C  and  allowed  to  reach  equilibrium.  Finally  the 
optics,  a  narrow  view  telescope,  was  attached  to  the  main  assemble.  The  narrow  view 
telescope  is  a  Cassegrainian  telescope  with  a  focal  length  of  12  cm.  It  has  a  narrow  field 
of  view  of  only  5  milliradians  or  0.286  degrees.  Cabling  was  connected  from  the 
spectrometer  to  the  computer.  This  completes  the  set  up  of  the  spectrometer. 

3.5.2  Spectrometer  Calibration 

Spectrum  calibration  was  conducted  following  the  procedures  in  the  Bomem  MR 
Series  Documentation  Set  (1997).  A  blackbody  was  place  in  front  of  the  spectrometer’s 
collimator  ensuring  it  filled  the  equipment’s  field  of  view.  This  essentially  eliminated  any 
stray  radiance.  The  blackbody  was  set  to  four  different  temperatures:  OC,  IOC,  30C  and 
50C,  and  allowed  to  stabilize  at  each  temperature.  Once  the  black  body  stabilized,  at  each 
temperature,  a  radiometric  reference  was  obtained.  After  the  radiometric  references  were 
obtained,  a  calibration  file  was  calculated  using  the  Compute  Radiance  Correction 
function  in  the  Research  Acquire  software.  This  completes  the  calibration  requirements. 

3.5.3  Acquiring  Atmospheric  Radiometric  Data 

After  the  setup  and  calibration  of  the  equipment,  acquiring  radiometric  data  is 
obtained  easily.  At  a  selected  time,  an  interferogram  is  obtained  using  the  spectrometer. 

A  interferogram  measurement  was  comprised  the  average  of  50  scans.  After  the 
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interferogram  is  obtained,  no  further  use  of  the  spectrometer  is  needed.  The  Research 
Acquire  software  constructs  radiance  file  from  the  interferogram,  then  the  calibrated 
reference  file  is  use  to  calibrate  the  radiance  file.  At  this  point,  calibrated  radiance 
measurements  can  be  taken  from  the  output  files  or  from  the  display  screen. 

3.6  Computing  Atmospheric  Temperature  Profiles 

The  key  to  this  method  is  to  determine  AB  (z),  the  change  in  brightness 
temperatures  as  a  function  of  height.  Once  AB  (z)  is  determined  it  is  used  to  calculate  the 
actual  brightness  temperatures  by  following  the  previous  assumption  that  the  actual 
brightness  temperatures  are  a  summation  of  the  brightness  temperature  used  in  PLEXUS 
and  a  change  in  the  brightness  temperature,  B(z)actuai  =  B(z)piexus  +  AB(z),  where  Bactuai  is 
the  physical  brightness  temperature  and  Bpiexus  is  the  brightness  temperature  used  by 
PLEXUS  software.  Temperature  is  calculated  from  Bactuai  using  Planck’s  Law. 

The  first  step  is  to  calculate  the  coefficients  for  AB(z).  This  is  accomplished  by 

using/,  -I2  =  J  W(z)AB(z)dz,  see  Equation  3.8.  W(z)  is  calculated  using  PLEXUS  derived 

kernel  functions.  The  change  in  radiance,  (Z,  -Zj),  is  calculated  using  spectrometer 
radiance  measurements,  Z,  and  PLEXUS  radiance  calculations,  Zj. 

The  Mathcad  templates,  calc_profile_2.mcd  and  calc_profile_3.mcd,  see 
Appendix  2,  are  used  to  retrieve  the  new  temperature  profiles.  Template 
calc_profile_2.mcd  assumed  AB(z)=a+bz  and  template  calc_profile_2.mcd  assumed 
AB(z)=a+bz+cz^. 
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IV.  Results 


4. 1  Introduction 

This  chapter  lists  the  results  observed  during  the  research  and  the  analysis  on 
success  of  the  physical  retrieval  method  employed.  First,  the  results  of  spectral  range 
determination  are  presented.  Next,  the  results  from  the  collection  of  radiometric  data  are 
presented.  Finally,  the  results  of  the  method  used  to  recover  atmospheric  temperature 
profiles  are  discussed. 

4.2  Spectral  Range  Determination 

Table  4.1  lists  the  wavenumbers  meeting  the  selection  criteria  presented  in 
Chapter  3.  First,  radiance  values  for  frequencies  influenced  by  highly  variable 
atmospheric  constituents  were  eliminated.  This  disqualified  a  significant  portion  of  the 
initial  650  cm’’  to  3000cm’’  region,  hence  wavenumbers  outside  of  the  2400  cm’’  and 
2800  cm’’  region  were  disqualified.  Secondly,  the  maximum  change  between  the  July  and 
January  transmittance  values  at  all  levels  was  required  to  be  less  than  five  percent.  Next, 
there  was  to  be  a  significant  change  between  the  January  and  July  radiance  values.  Figure 

3.3  was  used  to  subjectively  eliminate  the  region  between  2400  cm’’  and  2500  cm’’  . 
Finally,  kernel  functions  not  equal  to  zero  above  50,000  meters  were  eliminated.  This 
was  done  to  decrease  the  require  integration  of  the  radiative  transfer  equation  from  0  km 
-  300  km  to  0  km  -  50  km. 

Of  frequencies  listed  in  Table  4.1,  a  final  set  of  wavenumbers,  2566  cm’’  to  2570 
cm’’,  was  chosen  because  they  are  the  largest  set  of  continuous  wavenumbers  meeting  the 
criteria.  The  larger  set  allows  the  use  of  five  independent  pieces  of  information,  thereby 


28 


allowing  a  choice  in  the  assumed  change  in  brightness  temperature  function,  AB(z),  to 
have  up  to  five  coefficients. 


Table  4.1  Most  promising  frequencies. 


Wavenumber  (cm-1) 

Ratio  of 
January  and 
July 

Transmissivity 

values 

Height  where 
transmissivity 
change  with 
Height  =  0 
(km) 

Change  in 
Radiance  with 
height  =  0 
(km) 

Change  in 
Radiance 

Between  July  and 

2 

January  (W/cm 
sr  cm"') 

2534 

0.96038 

50 

50 

4.1655E-09 

2535 

0.97145 

50 

50 

3.3583E-09 

2536 

0.97491 

50 

50 

3.2157E-09 

2537 

0.97728 

50 

50 

3.0976E-09 

2538 

0.97460 

50 

50 

3.2201E-09 

2539 

0.96863 

50 

50 

3.7177E-09 

2540 

0.97797 

50 

50 

3.0406E-09 

2541 

0.97816 

50 

50 

3.0577E-09 

2542 

0.97875 

50 

50 

3.0566E-09 

2566 

0.97699 

50 

50 

2.5064E-09 

2567 

0.97666 

50 

50 

3.2534E-09 

2568 

0.97922 

50 

50 

2.7339E-09 

2569 

0.97636 

50 

50 

2.9552E-09 

2570 

0.95690 

50 

50 

4.2404E-09 

2609 

0.96152 

50 

50 

2.5490E-09 

2610 

0.96925 

50 

50 

2.0221E-09 

2611 

0.97059 

50 

50 

1.9233E-09 

2612 

0.96647 

50 

50 

2.1797E-09 

4.3  Collecting  Atmospheric  Radiance  Measurements 

Table  4.2  is  a  list  of  the  radiance  measurement  observed  by  the  employed 
spectrometer.  The  dates  were  chosen  for  clear,  cloudless,  night  conditions.  The  time  was 
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Radiance  (W/cm2  sr  cm-i)  Radiance  (W/cm2  sr  cm-i) 


chosen  to  correspond  to  radiosonde  data.  The  measurements  were  taken  on  two  dates,  1 
December  1999  at  OOOOZ  and  8  December  1999  at  OOOOZ. 


Figure  4.1  Calibrated  spectrometer  output  for  1  December  1999. 


Figure  4.2  Calibrated  spectrometer  output  for  8  December  1999. 


Table  4.2  Atmospheric  radiance  measurements  collected  by  the  spectrometer. 


Wavenumber  (cm  ’) 

Measured  Radiance  Values  at 
OOZ  1  December  1999 

(W/cm'^  sr  cm-1) 

Measured  Radiance  Values 
at  OOZ  8  December  1999 

(W/cm'^  sr  cm-1) 

2566 

1.300  X  10'’' 

1.721  X  lO'** 

2567 

1.384  X  10'" 

1.486  X  10'" 

2568 

1.491  X  10" 

1.712  X  10’’ 

2569 

1.199  X  10'" 

1.621  X  10'” 

2570 

1.601  X  lO  ** 

1.414X  10*’ 

4.4  Computing  Atmospheric  Temperature  Profiles 
4.4.1  Test  Case  I:  1  December  1999 

Figure  4.3  and  Figure  4.4  display  the  recovered  temperature  profile,  the 
radiosonde  temperature  profile  and  the  temperature  profile  used  in  calculating  the  kernel 
functions.  The  spectrometer  temperatures  in  Figure  4.3  were  calculated  by  assuming  the 
form  of  the  change  in  brightness  temperature,  AB(z),  was  AB(z)  =  a  +  bz.  The 
spectrometer  temperatures  in  Figure  4.4  were  calculated  by  assuming  the  form  of  the 
change  in  brightness  temperature,  AB(z),  was  AB(z)  =  a  +  bz  +  cz^  .  The  retrieval  of  a 
temperature  profile  was  successful,  although  the  recovered  temperature  profile  was 
significantly  different  than  the  radiosonde  derived  temperature  profile. 
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Temperature  (K)  Temperature  (K) 


Height  (m) 


Figure  4.3  Temperature  Plots  for  1  December  1999. 

.  Radiosonde  Temperatures _ PLEXUS  Temperatures 

- Spectrometer  Temperatures  using  AB(z)  =  a  +  bz. 


Height  (m) 


Figure  4.4  Temperature  Plots  for  1  December  1999. 

_ Radiosonde  Temperatures _ PLEXUS  Temperatures 

- Spectrometer  Temperatures,  using  AB(z)  =  a  +  bz  +cz  . 
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Brightness  Temeratures  (W/cnf  cm’'') 


Figure  4.5  displays  AB(z)  using  radiosonde  data  and  a  best-fit  function  for  that 
data.  Figure  4.6  displays  plots  of  AB(z)  using  radiosonde  data ,  AB(z)  =  a  +  bz  and  AB(z) 
=  a  +  bz  +  cz^. 


Figure  4.5  Plot  of  AB  using  1  December  1999  radiosonde  data. 


Height  (m) 


Figure  4.6  Plots  of  AB(z)  for  1  December  1999. 

AB(z)  from  radiosonde  data  . AB(z)  =  a  +  bz  - AB(z)  =  a  +  bz  4-  cz^. 
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4.4.2  Test  Case  II:  8  December  1999 


Figure  4.7  and  Figure  4.8  display  the  recovered  temperature  profile,  the 
radiosonde  temperature  profile  and  the  temperature  profile  use  in  calculating  the  kernel 
functions.  The  spectrometer  temperatures  in  Figure  4.7  were  calculated  by  assuming  the 
form  of  the  change  in  brightness  temperature,  AB(z),  was  AB(z)  =  a  +  bz.  The 
spectrometer  temperatures  in  Figure  4.8  were  calculated  by  assuming  the  form  of  the 
change  in  brightness  temperature,  AB(z),  was  AB(z)  =  a  +  bz  +  cz^.  Similar  to  test  case  I, 
the  retrieval  of  a  temperature  profile  was  successful,  although  the  recovery  of  the 
radiosonde  derived  temperature  profile  was  unsuccessful.  Figure  4.9  displays  AB  using 
radiosonde  data  and  a  best- fit  function  for  that  data.  Figure  4.10  displays  plots  of  AB(z) 
using  radiosonde  data,  AB(z)  =  a  bz  and  AB(z)  =  a  +  bz  +  cz^ 


Height  (m) 


Figure  4.7  Temperature  Plots  for  8  December  1999. 

. . . .  Radiosonde  Temperatures _ PLEXUS  Temperatures 

- Spectrometer  Temperatures  using  AB(z)  =  a  +  bz. 
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Brightness  Temperature  (\A//ci^  cm'^) 


Height  (m) 


Figure  4.8  Temperature  Plots  for  1  December  1999. 

....  Radiosonde  Temperatures _ PLEXUS  Temperatures 

- Spectrometer  Temperatures,  using  AB(z)  =  a  +  bz  +cz  . 


5e-08 
O 

-5e-08 
-1e-07 
-1 .5e-07 
-2e-07 
-2.5e-07 

0  2500  5000  7500  10000  12500 

Height  (m) 


Figure  4.9  Plot  of  AB  using  8  December  1999  radiosonde  data. 


5.0E-07 


Height  (m) 

Figure  4.10  Plots  of  AB(z)  for  8  December  1999. 

_ AB(z)  from  radiosonde  data  . AB(z)  =  a  +  bz  - AB(z)  =  a  +  bz  +  cz^. 

4.4.3  Summary 

The  fitting  of  functions  to  the  actual  change  in  brightness  temperatures  revealed 
the  better  form  of  AB(z )  to  be  a  polynomial  of  higher  order  than  used  in  this  research. 
Computer  precision  in  calculating  the  kernel  matrix  and  its  inverse  limited  the  choice  of 
assumed  forms  of  AB(z)  to  no  larger  than  a  second-degree  polynomial.  A  qualitative 
analysis  of  the  recovered  temperatures  profiles  and  the  assumed  forms  of  AB(z  )  reveals 
the  assumed  form  of  AB(z )  needs  to  be  a  higher  order  polynomial.  In  both  cases,  the 
AB(z )  functions  using  the  higher  order  polynomial  produced  a  temperature  profile, 
which  is  closer  to  the  radiosonde  derived  temperature  profile  for  the  lower  atmosphere. 
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IV.  Conclusions  and  Recommendations 


The  method  employed  in  this  research  was  able  to  recover  an  atmospheric 
temperature  profile,  however  significant  differences  were  observe  when  the  recovered 
temperature  profile  was  compared  to  a  temperature  profile  recovered  by  a  radiosonde  for 
the  same  time  frame.  The  investigation  of  possible  choices  of  frequencies  to  use  in  the 
method  revealed  various  frequencies.  Due  to  the  fact  that  highly  variable  atmospheric 
constituents  affected  the  radiance  values  for  most  of  the  initial  region,  the  2500  cm"'  to 
2800  cm‘'  region  displayed  the  most  suitable  frequencies. 

Using  functions  fitted  to  model  data,  kernel  and  brightness  temperature  functions, 
spectral  radiances  were  calculated  using  Schwarzchild's  equation.  The  calculated  spectral 
radiances  did  not  agree  with  the  spectral  radiances  calculated  from  the  PLEXUS 
software.  Further  research  is  needed  to  explain  the  differences. 

The  method  researched  demonstrated  its  ability  to  recover  an  atmospheric 
temperature  profile  and  I  recommend  further  research  of  this  method.  A  greater  accuracy 
was  observed  when  the  assumed  form  of  AB(z),  change  in  brightness  temperature 
function,  used  higher  order  terms.  Computer  precision  limited  the  assumed  polynomial 
forms  of  AB(z)  to  second-degree  polynomials.  An  interesting  follow-on  research  project 
would  be  to  explore  the  effects  of  using  different  forms  for  the  brightness  temperature 
function,  kernel  functions,  and  AB(z)  functions  and  how  this  affects  the  accuracy  of  the 
recovered  temperature  profiles.  The  condition  numbers  of  the  matrices  are  large,  on  the 
order  of  10^,  and  indicate  the  matrices  may  be  ill  conditioned.  The  use  of  different  forms 
for  the  brightness  temperature  function,  kernel  functions  and  AB(z)  may  improve  the 
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accuracy  of  the  recovered  temperature  profile  and  reduce  the  condition  numbers  of  the 
matrices. 

In  the  calibration  process,  the  blackbody’s  minimum  temperature  was  zero 
degrees  Celsius.  This  affected  the  calibration  files  used  by  the  BOMEM  software  as 
minimum  atmospheric  temperatures  are  well  below  zero  degrees  Celsius.  I  recommend 
obtaining  a  blackbody  capable  of  lower  minimum  temperatures  to  use  in  the  calibration 
process.  In  verifying  the  recovered  temperatures,  the  radiosonde  and  spectrometer  sites 
were  approximately  thirty  miles  apart.  If  possible,  the  radiosonde  and  spectrometer  sites 
should  be  the  same  to  help  reduce  error. 
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Appendix  1:  Brightness  Temperature  and  Transmittance  Function  Data 
This  section  contains  the  best-fit  function  data  for  brightness  temperature  function 
as  well  as  the  transmittance  and  kernel  parameters  for  2566  cm  ’  to  2570cm'’. 


39 


Brightness  Temperature  Data 

B(z)=(a+cz+ez^+gz^+iz''+kz^)/(l+bz+dz^+fz^+hz'’+jz^)  is  the  function  chosen  to  represent 
the  brightness  temperatures  used  by  the  PLEXUS  software  to  calculate  transmittance  and 
radiance  values.  Table  A.l  .1  list  the  values  and  statistical  data  of  the  parameter  for  the 
equation.  The  adjusted  r2  value  for  the  fitted  equation  is  0.99999994. 

Table  A.  1.1  Numeric  summary  of  fitted  brightness  temperature  function  parameters. 


Farm 

Value 

Std_Error 

t-value 

P>|t| 

a 

3.90765E-07 

4.16314E-11 

9386.29893 

0.0000 

b 

-2.00020E-04 

4.34230E-06 

-46.06311 

0.0000 

c 

-1.17790E-10 

1.73953E-12 

-67.71151 

0.0000 

d 

5.95852E-08 

1.13240E-09 

52.61859 

0.0000 

e 

1.58601E-14 

3.42736E-16 

46.27508 

0.0000 

f 

-8.32550E-12 

3.61496E-13 

-23.03055 

0.0000 

g 

-1.13980E-18 

3.881 07E-20 

-29.36751 

0.0000 

h 

7.99827E-16 

7.85057E-17 

10.18815 

0.0000 

i 

4.30596E-23 

2.64977E-24 

16.25029 

0.0000 

j 

-2.74370E-20 

5.37737E-21 

-5.10239 

0.0005 

k 

-6.78200E-28 

7.27836E-29 

-9.31799 

0.0000 
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Brightness  Temperature  (m  W/m  2)  Brightness  Temperature  (m  w/m  2) 


y=(a+cx+ex2+gx^+ix^+kx^)/(1  +bx+d)^+fx2+hx^+jx?) 


y=(a+cx+ex  2+gx^+ix'^+kx^)/(1 +bx+dx  2+fx^+hx^+jx^) 


Figure  A.  1 .2  Residual  plot  of  fitted  brightness  temperature  function. 
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2566  cm  *  Data 

't2566(z)=(a+cz)/(l+bz+dz^)  is  the  function  chosen  to  represent  the  transmittance 
values  calculated  by  the  PLEXUS  software.  Table  A.  1.2  lists  the  values  and  statistical 
data  of  the  parameter  for  the  equation.  The  adjusted  r2  value  for  the  fitted  equation  is 
0.9995. 

Table  A.  1.2  Numeric  summary  for  transmittance  and  kernel  parameters  for  2566  cm 


Parm 

Value 

Std_Error 

t-value 

P>|t| 

a 

0.999910333 

9.971 47e-05 

10027.70781 

0.0000 

b 

0.000308564 

8.31176e-07 

371.2381899 

0.0000 

c 

0.000262176 

7.32808e-07 

357.7684234 

0.0000 

d 

-1.1198e-10 

4.51076e-13 

-248.259375 

0.0000 
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Transmittance 


y=(a+cx)/(1  +bx+dx2) 
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Figure  A.  1 .4  Residual  plot  of  fitted  transmittance  function  for  2566  cm 


2567  cm  ’  Data 

T2567(z)=(a+cz)/(l+bz+dz^)  is  the  function  chosen  to  represent  the  transmittance 
values  calculated  by  the  PLEXUS  software.  Table  A.  1.2  lists  the  values  and  statistical 
data  of  the  parameter  for  the  equation.  The  adjusted  r2  value  for  the  fitted  equation  is 
0.9997224799. 

Table  A.  1.3  Numeric  summary  for  transmittance  and  kernel  parameters  for  2567  cm"’. 


Farm 

Value 

Std_Error 

t-value 

P>|t| 

a 

0.998440828 

0.00011731 

8511.142627 

0.0000 

b 

0.000206724 

5.15083e-07 

401.3422029 

0.0000 

c 

0.000151993 

4.19289e-07 

362.500531 

0.0000 

d 

-2.1125e-10 

6.49634e-13 

-325.190357 

0.0000 
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Transmittance 


y=(a+cx)/(1+bx+clx2) 


Figure  A.  1 .5  Plot  of  transmittance  values  and  fitted  function  for  2567  cm’'. 


y=(a+cx)/(1+bx+clx2) 


Figure  A.  1 .6  Residual  plot  of  fitted  transmittance  function  for  2567  cm’' . 
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2568  cm'^  Data 

T^2568(z)^(a+cz)/(l+bz+dz^)  is  the  function  chosen  to  represent  the  transmittance 
values  calculated  by  the  PLEXUS  software.  Table  A.  1.2  lists  the  values  and  statistical 
data  of  the  parameter  for  the  equation.  The  adjusted  r2  value  for  the  fitted  equation  is 
0.9998. 

Table  A.  1.4  Numeric  summary  for  transmittance  and  kernel  parameters  for  2568  cm'’. 


Farm 

Value 

Std_Error 

t-value 

P>|t| 

a 

0.999895747 

7.8544e-05 

12730.3925 

0.0000 

b 

0.000237587 

4.45246e-07 

533.6068558 

0.0000 

c 

0.00018786 

3.76929e-07 

498.3963266 

0.0000 

d 

-1.6576e-10 

3.95142e-13 

-419.50688 

0.0000 
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Figure  A.  1 .7  Plot  of  transmittance  values  and  fitted  function  for  2568  cm"' . 
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Figure  A.  1 .8  Residual  plot  of  fitted  transmittance  function  for  2568  cm'’ . 


47 


2569  cm  Data 


T2569(z)=(a+cz)/(l+bz+dz^)  is  the  function  chosen  to  represent  the  transmittance 
values  calculated  by  the  PLEXUS  software.  Table  A.  1.2  lists  the  values  and  statistical 
data  of  the  parameter  for  the  equation.  The  adjusted  r2  value  for  the  fitted  equation  is 
0.9998292345. 


Table  A.  1.5  Numeric  summary  for  transmittance  and  kernel  parameters  for  2569  cm’’. 


Parm 

Value 

Std_Error 

t-value 

P>|t| 

a 

1.000358467 

8.185e-05 

12221.85763 

0.0000 

b 

0.000243899 

4.48959e-07 

543.2543509 

0.0000 

c 

0.000190951 

3.7701 1e-07 

506.4856435 

0.0000 

d 

-1.7589e-10 

4.10406e-13 

-428.583981 

0.0000 

48 


Transmittance  Transmittance 


Figure  A.  1.9 


Plot  of  transmittance  values  and  fitted  function  for  2569  cm’'. 
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Figure  A.  1 . 1 0  Residual  plot  of  fitted  transmittance  function  for  2569  cm’' . 
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2570  cm  Data 


't257o(z)=(a+cz)/(l+bz+dz^)  is  the  function  chosen  to  represent  the  transmittance 
values  calculated  by  the  PLEXUS  software.  Table  A.  1.2  lists  the  values  and  statistical 
data  of  the  parameter  for  the  equation.  The  adjusted  r2  value  for  the  fitted  equation  is 
0.9998559591 . 


Table  A.  1.6  Numeric  summary  for  transmittance  and  kernel  parameters  for  2570  cm"’. 


Parm 

Value 

Std_Error 

t-value 

P>|t| 

a 

1.001225275 

8.88782e-05 

11265.13557 

0.0000 

b 

0.000289502 

4.5391 3e-07 

637.7904422 

0.0000 

c 

0.000219657 

3.69277e-07 

594.8294418 

0.0000 

d 

-2.081 4e-10 

4.35863e-13 

-477.529018 

0.0000 
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APPENDIX  2:  IDL  Programs  and  Mathcad  Templates 
Two  IDL  programs  and  two  Mathcad  templates  are  provided. 

The  two  IDL  programs  are  Plexusbatch.pro  and  ADI.pro.  Plexusbatch.pro 
increments  the  pathlengths  in  the  PLEXUS  input  file  then  call  the  PLEXUS  program  to 
compute  radiance  and  transmittance  value  for  that  pathlenght.  ADI.pro  manipulates 
PLEXUS  outputs.  It  takes  200  individual  transmittance  file  produced  by  PLEXUS  and 
put  them  into  one  files  in  a  column  format  for  easy  importation  into  Microsoft  Excel. 

The  two  Mathcad  templates  are  Recover_Profile_2.MCD  and 
Recover_Profile_3.MCD.  Recover_Profile_2.MCD  assumes  AB(z)  =  a  +  bz,  solves  for 
AB(z)  and  calculate  the  a  temperature  profile.  Recover_Profile_3.MCD  assumes  AB(z)  = 
a  +  bz  +  cz^,  solves  for  AB(z),  and  calculate  the  a  temperature  profile. 
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pro  plexusbatch 


Plexusbatch.pro 


;  This  program  increments  the  pathlengths  in  the  PLEXUS  input  file  then  call  the 
PLEXUS  program  to  compute  radiance  and  transmittance  value  for  that  pathlenght. 


;  closes  all  devices  left  open 
for  j= 1,200  do  begin  ;  loopl 
close, /all 

;  allows  selection  of  a  input  and  output  files 


fin- c:\plexus\data_out\adi\height.txt' 
fh- c  :\plexus\adi\adi  .cpd' 


;  defines  a  string  variable 


;  returns  the  number  of  rows  in  the  file 

while  not  (EOF(2))  do  begin 
readf,2,s 

n=n+l 

endwhile 

;total=n-l 

close,2 

;  read  pathlength  line  of  input  file 

openr,2,fn 
lines=strarr(n) 
readf,2, lines 
close,2 

var=j*0.010 

var=strcompress(string(var)) 

lines(29)="Final_Altitude="+var+"  ;[H2ALT-F,M,S,U]  (KM)" 
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file_name='c:\plexus\adi\adi.cpd' 
file_name2- c:\plexus\data_out\adi\height.txt' 
close,  /all 

;  Writes  new  pathlength  line  to  PLEXUS  input  file; 

openu, 2, file_name;, /append 
openu,3,file_name2, /append 
for  i=01,n-l  do  begin 
;print,  i 

printf,2,lines(i) 


;  Writes  pathlength  value  to  file; 

Endfor 

printf,3  ,j ,  var*  100.0 
print,  'wrote  the  file' 

close,/all 

spawn,  'plexus.exe  c:\plexus\adi\adi.cpd' 
endfor 
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ADI.PRO 


pro  adi 

;  This  program  takes  200  individual  transmittance  file  produced  by  PLEXUS  and  put 
them  into  one  files  in  a  column  format. 

;  closes  all  devices  left  open 

close, /all 

;allows  selection  of  a  specific  file 


fh='c:\plexus\adi\excel2  .out' 


output=fltarr(20,235 1 ) 
for  i=0,19  do  begin 

num=strcompress(sindgen(1000),/remove_all) 
num(0: 9)='0'+num(0 : 9) 
num(0: 99)='0'+num(0: 99) 
var=num(i+431) 
print,  var 


fin- E:\l  5jul99\adi'+var+'.tm' 

print,  fin 

openr,2,fm 

data=fltarr(2,235 1 ) 

readf,2,data 

close,2 

output(i,*)=data(  1  ,*) 
endfor 

openw,3,fn,  width=1300 
format='(20(E  1 2 . 5 , 1  x))' 
printf,3  ,output,format=format 
close,2 
close,/all 
end 
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Recover  Profile  2.MCD 


ORIGII^I 

This  mathcad  template  calculates  a  temperature  profile  assuming  the  change  in  the 
brightness  temperature  function  is  A  B(z)  :=  +  b  . 

Assuming  B  is  linear  in  form 
AB(z)  ;=  (a^  +  bjj-z) 

Transmittance  and  kernel  functions  for  2567  cm~l 
a  256?:=  0.99844082; 
b  2557;= 0.00020672' 
c  2567-0.00015199; 
d2567'=-2-l  12510-"' 

,  ,  ^  2567+  2567^ 

'  2561<^)  == - J 

1+1^  2567^+ ‘1 2567^ 

,  ,  ■  ^  2567+  2567^  2567^^  +  ^  25671’  2567+  2'a  2567*1  2567^) 

W  256'^ z)  --2 - - - - - 

(l  +  l’2567^+^2567M 


Transmittance  and  kernel  functions  for  2568  cm'l 

a  2568  ■=0.99989574’ 
b  2568-0.00023758' 
c  2568- 0.0001 878( 
d2568’=-l-657610-'« 
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^  256^2)  ==- 


**  2568+  2568'^ 


W256^") 


'  +  *’  2568'^+ ‘^2568’^ 

2568+  2568*^  2568’^^  +  ^ 2568*^  2568+  2568'^  2568’^, 

(l  +  *’ 2568'^+ ‘^  2568’^^) 


Solving  for  kernel  functions  coefficients 

•55000 

AB(z)  W  2567(z)dz  simplify  .2007062487288491 151d  j,  -  828.130935443368960019 ^ 

Jo 


•55000 

AB(z)  W  256^z)dz  simplify  -.1645491 555083 136263Sj,-  645.601385564716100019 

Jo 


W 


Kernel  Function  Coeff. 

.2007062487288491151 1828.1309354433689600; 
.  1 64549 1 5 55083 1 362639645 .60 1 3  8 55647 1 6 1 000 ; 


Radiance  values  for  2567  cm-1  and  2568  cm-1,  respectively,  from  PLEXUS  output. 


'  plexus 


7.48100010'^ 

6.19910010'^ 


Radiance  values  for  2567  cm-1  and  2568  cm-1,  respectively,  obtained  using 
spectrometer. 


'  spect 


1. 38610** 
1. 491' 10** 


•weight 
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Calculating  (L^  -  L2) 

A  L  j  .  -  L  spgcj  —  L  plexus 

Radiance  difference  for  2567  cm-l  and  2568  cm-l,  respectively. 

-1.0776810'’ 

ALj  = 

6.89320000000002  O'*’ 


Solving  for  the  change  in  brightness  temperature  function.  DBtz).  coefficients 


-I 


coeff  1  :=W  -ALj 


coeff  j  = 


1.8927021883187-20“ 


4.7172969811734-30 


11 


a b -coeff. 


bb- coeff. 


ab  =  -l-8927021883187-20 


r7 


bb  =  4.717296981173430 


rll 


AB  j(z)  ;=ab-l-bb'Z 

Brightness  temperature  function  used  to  calculate  initial  kernel  functions 
a  b  :=  3.9076510"’^ 
bb  --2.OOO2IO"’ 
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;=-1.177910‘'° 
d^,  ;=5.9585210'* 
1.586 10' 

f^,  ;=-8.325510‘‘^ 

1 R 

gb:=-1.1398ld'‘'* 
;=7.9982710‘'® 
i^,  ;=4.3059610‘^^ 
jb  :=- 2.7437 10'^® 
;=- 6.782 10'^* 


®plexus(^)  ■  2  3  4  5 

l  +  bt,  +  b|jZ+d|jZ  d-f^z  +h^,z  +i^,z 


i;=0..  55 

Calculating  Bfz)  +  AB(z) 

Bl  ■=(Bplexus(l  >»«»  +  *B,(i  lOOO)).[em-'.r5E] 

\  cm  / 

Converting  new  brightness  temperatures  into  a  temperature  profile 

T 

^  spectrometer 
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T 


spectrometer 


1 

1 

274.82 

2 

274.21 

3 

271.75 

4 

268.97 

5 

267.3 

6 

267.53 

7 

269.52 

8 

272.57 

9 

276.03 

10 

279.51 

11 

282.8 

12 

285.86 

13 

288.68 

14 

291.26 

15 

293.65 

16 

295.87 

K 


‘290.233 

'274.8 

286.949 

266.2 

281.901 

262.2 

275.579 

261,0 

268.449 

257.7 

260.924 

K  T 

radiosonde 

252.6 

253.355 

245.9 

246.023 

240.0 

239.15 

233.7 

232.909 

226.3 

227.433 

_218.1_ 

Calculated  Temperatures 
using  spectrometer  method 


Plexus  inputs  Actual  Temperatures 

from  radiosonde 


h  10 

Calculate  difference  between  calulated  temperature  profile  and  radiosonde  temperature  profile. 


AT  :=T 

spectrometer 


T 

radiosonde. 


1 

1 

0.019066430679402 

2 

8.0123603635991 

3 

9.55369440331771 

4 

7.96902201582952 

5 

9.60055267729683 

6 

14.9324461347775 

7 

23.6155693849379 

8 

32.5668106084541 

9 

42.3318372053101 

10 

53.207175648985 
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Recover  Profile  3.MCD 


ORIGIT^l 

This  mathcad  template  calculates  a  temperature  profile  assuming  the  change  in  the 
brightness  temperature  function  isAB(z)  ;=a^,-t-bj,  z+c  . 

Assuming  B  is  linear  in  form 

AB(z)  ;=aj,-+-bj,z+C|jZ^. 

Transmittance  and  kernel  functions  for  2567  cm~^ 
a  256?;=  0.99844082: 
b  2567- 0.00020672' 
c  2567-0.00015199; 
d2567--2.1 12510"'“ 

,  ,  2 2567+ ^2567^ 

T  256^z)  := - ^ 

1  +  ^2567^+*^  2567^ 

•  2567+  2567^*  2567^^+  ^  2567^’  2567+  ^  ^  2567*^  2567^) 

W  256^2)  - 

(’  +  *’2567^+ *^2567^^) 

Transmittance  and  kernel  functions  for  2568  cm~l 

a  2568-0.99989574- 
b  2568-0  00023758' 
c  •“  0.0001 878( 

d2568--l-657610-'“ 
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^  256^2)  == 


^  2568+  2568'^ 


W  256^") 


^  +  2568^+ ‘*2568'^ 

,_•  (•‘^  2568+  ^  2568’^  2568'^^  +  ^  2568^  2568+  2568'^  2568'^; 

2\^ 

*  +  2568^+  ^  2568^  j 


Transmissivity  and  Kernel  functions  for  2569  cm'^ 


2. 2569 ^  *00035846 


b  2569 ‘=0*00024389' 


c  2569 '=0.000 19095 


d2569-l-758910-'‘» 


,  ,  ^2569+^^2569'^ 

T  256c(z)  := - ^ 

*+^25692+*^  2569^ 


f  2  ^ 

■  H  2569+  ^  2569^^  2569^  +  ^  2569'’’  2569+  ^  ^  2569^^  2569^^ 

W  256s(z)  :=^ - - - - 

('  +  ’’2569^+‘^2569M 


Solving  for  kernel  functions  coefficients 

•55000 

A  B(z)-W  2567  (z)dz  simplify  -^  106132.17125040000000  -c  ^  - -2007062487288491 151 1  -36-  828.13093544336896000  -b  ^ 

Jo 

•55000 

AB(z)  W256^z)dz simplify  ^-197940.9267828000006^-  .1645491555083136263%- 645.6013855647161005% 

Jo 

•55000 

AB(z)  W 256^z)dz simplify  “^-18185.15202290000006%-  .1717782573165399743%-  662.1342997687510005% 

Jo 
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Weighting  Function  Coefficients 


w 


.2007062487288491 1 5 1 1828. 1 30935443368960001 06 1 32.1 7 1 2504000000 : 
.  1 64549 15550831 362639645.60 1 3855647 1 6 1 00001 97940.9267828000000 : 
.171 7782573 1 653997457562. 1 3429976875 1 000001 8185.1 5202290000000 ! 


Radiance  values  for  2567  cm-1, 2568  cm-l  and  2569  cm-l,  respectively,  from  PLEXUS 
output. 


5.52430010'^ 


"plexus  *“ 


4.54500010'^ 


4.79820010'^ 


Radiance  values  for  2567  cm-l,  2568  cm-l  and  2569  cm-l,  respectively,  obtained  using 
spectrometer. 


1. 38410"® 


'  spect 


L49MCr 


•weight 


1.19910"® 


Calculate  change  in  Radiance 
A  L  j  -  L  gpg(,(  —  L  p]exus 

Radiance  difference  for  2567  cm-l,  2568  cm-l  and  2569  cm-l,  respectively. 
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9.32602  Id" 


AL  j  = 


1.14534310'^ 


8.0670710'^ 


Solving  for  the  change  in  brightness  temperature  function,  DB(z).  coefficients 

coeffj  :=W'’  aLj 


-1.896968963650700'’ 


coeff  I  = 


3.766163063027040'“ 


-  2.3003665245562  0'’"^ 


a  1^1  :=  coeff  | 
bbi  :=  coeffj^ 

Cbl  ==  coeffj^ 

a^i  =-1.896968963650700'’ 
bbi  =  3.766163063027040'“ 
Cbl  =-2.30036652455620'“ 


ABi(z)  :=abi+bbi  z-|-Cbi  z’ 


Brightness  temperature  function  used  to  calculate  initial  kernel  functions 

ab  ;=  3.90765 10'’ 
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:=- 2.000210“* 
ci,:=-1177910‘‘® 
dj,  :=  5.9585210'* 
ej,  :=  1.586 10' 
fj,  :=-8.325510''^ 

1  R 

gj,:=-1.139810'‘'* 
h^,  :=  7.99827 10" 
ij,  :=  4.30596 10'^* 
jb:=-2.743710'^‘’ 
kb:=- 6.782 10'^* 

®plexus^^)  ■  2  3  4  5 

l  +  bj,  +  bjj.z4-d|j.z  +fbz  +lib'^  +jb'^ 


i:=0.,  55 

Calculating  B(z)  +  AB(z) 

Bl.^  ==(Bp|exus(ilOOO)  +  ABi(ilOOO))L-'-^ 

\  cm 

Converting  new  brightness  temperatures  into  a  temperature  profile 

T 

spectrometer. 


'Bi 
\  1+1/ 
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spectrometer 


1 

1 

274.82 

2 

274.21 

3 

271.75 

4 

268.97 

5 

267.3 

6 

267.53 

7 

269.52 

8 

272.57 

9 

276.03 

10 

279.51. 

11 

282.8 

12 

285.86 

13 

288.68 

14 

291.26 

15 

293.65 

16 

295.87 

K 


plexus 


*290-233 

*274.8' 

286.949 

266.2 

281.901 

262.2 

275.579 

261.0 

268.449 

257.7 

260.924 

radiosonde  ‘ 

252.6 

253.355 

245.9 

246.023 

240.0 

239.15 

233.7 

232.909 

226.3 

227.433 

.218.1. 

K 


Calculated  Temperatures 
using  spectrometer  method 


Plexus  inputs  Actual  Temperatures 

from  radiosonde 


h  :=l..  10 

Calculate  difference  between  calculated  temperature  profile  and  radiosonde  temperature  profile. 


spectrometer  radiosonde 


1 

1 

0.019066430679402 

2 

8.0123603635991 

3 

9.55369440331771 

4 

7.96902201582952 

5 

9.60055267729683 

6 

14.9324461347775 

7 

23.6155693849379 

8 

32.5668106084541 

9 

42.3318372053101 

10 

53.207175648985 
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APPENDIX  3:  Radiosonde  Data 


This  appendix  contains  radiosonde  data  for  the  two  test  cases.  Radiosonde  data  is  from 
the  Wilmington  OH  radiosonde  site,  which  is,  located  approximately  30  miles  southeast 
from  were  the  radiometric  data  was  collected. 

1  December  1999  QOOOZ  RAPE 


TTAA  51001  72426  99998  00058  01007  00309  Hill  Hill  92926  05956  35513  85581 
12110  34513  70080  10965  35537  50561  24759  35059  40719  37563  35570  30910  54357 
01086  25024  62556  36079  20163  57961  34076  15345  55173  34047  10604  57175  33029 
88239  63156  35067  77307  01087  40805  51515  10164  00021  10194///// 35529= 


UJXX  OOKW  010000  72426 

TTBB  51000  72426  00998  00058  1 1987  00960  22850  121 10  33830  13510  44822  08564 
55820  08375  66748  07988  77678  11362  88612  14563  99581  17964  11565  19960  22517 
22962  33462  29756  44400  37563  55279  58556  66239  63156  77198  57162  88173  58967 
99159  54972  11134  54176  22100  57175  31313  45202  82307  41414  856//= 


8  December  1999  QOOOZ  RAOB 


USXX  OOKW  080000  72426  TTAA  58001  72426  99986  02247  17005  00209  Hill  Hill 
92838  02668  Hill  85527  07493  23023  70107  00689  23539  50571  18179  25038  40734 
30723  25562  30932  45750  25583  25051  52562  26087  20193  57369  26575  15375  58373 
26551  88189  58368  26077  77290  25093  42509  51515  10158  10164  00014  10194  Hill 
24029= 

UJXX  OOKW  080000  72426  TTBB  58000  72426  00986  02247  1 1975  03458  22935 
00656  33925  02668  44917  04674  55860  07296  66850  07493  77837  06873  88739  04286 
99478  20977  11460  22966  22442  24736  33422  27316  44389  32321  55383  32157  66376 
33160  77353  37156  88323  42358  99285  48737  11276  50547  22250  52562  33189  58368 
44155  56775  55133  58974  66/// ///// 77123  61374  88120  61774  31313  45202  82302 
41414  00900  51515  10150  10158= 
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